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In this paper we propose some very promissing results in interval arithmetics which permit to build 
well-defined arithmetics including distributivity of multiplication and division according addition 
and substraction. Thus, it allows to build all algebraic operations and functions on intervals. This 
will avoid completely the wrapping effects and data dependance. Some simple applications for matrix 
eigenvalues calculations, inversion of symmetric matrices and finally optimization are exhibited in 
the object-oriented programming language python. 
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I. HISTORY 

The first mathematician who has used intervals was the famous Archimedes from Syracuse (287-212 
b.C). He has proposed a two-sides bounding ofTr: 3-|-^<7r<3-|-| using polygons and a systematic 
method to improve it. In the beg inning of the twentieth century, the American mathematician and physicist 
Wiener, published two papers^EI, and used intervals to give an interpretation to the position and the time 
of a system. More papers on the subject were writterP^^ only after Second World War. Nowadays, we 
consider R.E. Moore^li^ as the first mathematician who has proposed a framework for interval arithmetics 
and analysis. The interval arithmetic, or interval analysis has been introduced to compute very quickly 
range bounds (for example if a data is given up to an incertitude). Now interval arithmetic is a computing 
system which permits to perform error analysis by computing mathematic bounds. The extensions of the 
areas of applications are important: non linear problems, PDE, inverse problems. It finds a large place of 
applications in controllability, automatism, robotics, embedded systems, biomedical, haptic interfaces, form 
optimization, analysis of architecture plans, ... 

Interval calculations are used nowadays as a powerful tool for global optimization and set inversiorP^. Several 
groups have developped some softwares and libraries to perform those new apparoaches such as INTLAB^, 
INTOPT90 and GLOBSOL^,Numericc^. But our goal for this article, is not to replace the semantic 
approach of intervals, which has to be adapted to each problem by the engineer or the scientist, but to 
propose a new arithmetic of intervals, which allows to avoid the wrapping and data dependance effects. It 
yields to a better construction of inclusion functions. 

We expose in this paper the main results of a PhD thesis!' defended by one the author and some consecutive 
numerical applications. The plan is the following. In a first time we define a real Banach structure on 
the completion IM of the semi group of intervals IK, with a vector space structure. This permits to define 
the notion of differential function with values of IM and to use some important tools and the fixed point 
theorem. Next we extend the classical product to have a distributivity property. With this approach we 
obtain a notion of differential calculus and a natural linear algebra on the set of intervals. After that, we 
gives some examples in a python implementation and we end this article by giving some simple numerical 
applications : optimization of interval functions, interval matrix diagonalization, and inversion of symmetric 
matrices . 



II. AN ALGEBRAIC APPROACH TO THE SET OF INTERVALS. 

In this section we present the set of intervals as a normed vector space. We define also a four-dimensional 
associative algebra whose product gives the product of intervals in any cases. 

A. Minkowski operations 

An interval is a bounded non empty connected closed subset of M. Let IE be the set of intervals. The 
semantical arithmetic operations on intervals, called Minkowski operations, are defined such that the result 
of the corresponding operation on elements belonging to operand intervals belongs to the resulting interval. 
That is, if o denotes one of the semantical operations -I-, — , *, we have, if X and Y are bounded intervals of 

M, 

XoY = {xoy / xex, yev}, 

In many problems using interval arithmetic, that is the set IM with the Minkowski operations, there 
exists an informal transfers principle which permits, to associate with a real function / a function define 
on the set of intervals IM which coincides with / on the interval reduced to a point. But this transferred 
function is not unic^ue. For example, if we consider the real function f{x) = x'^ + x = x(x + 1), we associate 
naturally the functions /i : IM — > M given by fi{X) = X{X + 1) and hiX) ^ X'^ + X. These two 
functions do not coincide. Usually this problem is removed considering the most interesting transfers. But 
the qualitative "interesting" depends of the studied model and it is not given by a formal process. In this 
section, we determine a natural extension IM of IM provided with a vector space structure. The vectorial 
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substraction X \ y does not correspond to the semantical difference of intervals and the interval \X has 
no real interpretation. But these "negative" intervals have a computational role. If a problem conduce to a 
"negative" result, then this problem is "pervert" (see Lazare Carnot with his feeling on the natural negative 
number) . 

Let IM be the set of intervals. It is in one to one correspondence with the half plane of M^: 

Vi = {{a,b),a<b}. 

This set is closed for the addition and Vi is endowed with a regular semi-group structure. Let V2 be the 
half plane symmetric to Vi with respect to the first bisector A of equation y — x = 0. The substraction on 
IM, which is not the symmetric operation of +, corresponds to the following operation on Vi: 

(a, b) - (c, d) = (a, b) + sa° so(c, d), 

where sq is the symmetry with respect to 0, and sa with respect to A. The multiplication * is not globally 
defined. Consider the following subset of Vi : 

■Pi.i ={(a,6) eri,a>0,b> 0}, 
Vi.2 = {{a,b)Gri,a<0,b>0}, 
Vi,3 = {ia,b) £ Vi,a < 0,b < 0}. 

We have the following cases: 

1) If (a, &), (c, d) € Pli the product is written (a, b) * (c, d) = {ac, bd). 

The vectors ei — (1, 1) and 62 — (0, 1) generate T^i^i that is any {x, y) in 7^i,i, can be decomposed as 

(x, y) = xei + {y — x)e2, with x > and y ~ x > 0. 
The multiplication corresponds in this case to the following associative commutative algebra: 

eiei = ei, 

6162 = 6261 — 6262 = 62. 

2) Assume that (a, b) e T^i.i and (c, d) e 7^i,2 so c < and d > 0. Thus we obtain (a, b) * (c, d) — {be, bd) 
and this product does not depend of a. Then we obtain the same result for any a < b. The product 
(a, b) * (c, d) — {be, bd) corresponds to 

eiei = 6261 = ei 

6162 = 6262 = 62 

This algebra is not commutative and it is different from the previous. 

3) If (a, b) e T'l.i and (c, d) £ Vi^s then a > 0, 6 > and c < 0, d < and we have (a, b) * (c, d) — {be, ad). 
Lot 61 = (1, 1), 62 = (0, 1). This product corresponds to the following associative algebra: 




-62. 

This algebra is not associative because (6261)61 7^ 62(6161). We have similar results for the cases 
(^1,2,^1,2), (^1,2, 7^1,3) and (751,3, ^1,3). 
An objective of this paper is to present an associative algebra which contains all these results. 



B. The real vector space IR 

We recall briefly the construction proposed by Marko\Elto define a structure of abelian group. As (IM, +) 
is a commutative and regular semi-group, the quotient set, denoted by IM, associated with the equivalence 
relations: 



{x, y) r^ {z, t) <=^ x + t = y + . 
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for all x, y, z,t IE IM, is provided with a structure of abelian group for the natural addition: 



{x,y) + {z,t) = {x + z,y + t) 



where (x, y) is the equivalence class of {x, y). We denote by \(.t, y) the opposite of {x, y) . We have \(a;, y) = 
{y,x). li X = [a, a], a G M, then (x, 0) = (0, — x) where —x ~ [—a,— a], and \(a;,0) = (0,a;). In this case, 
we identify x — [a, a] with a and we denote always by M the subset of intervals of type [a, a]. Naturally, the 
group IE is isomorphic to the additive group M^ by the isomorphism (([a, 5], [c, d]) — >■ (a — c, & — d). We find 
the notion of generalized interval. 

Proposition 1 Let X = {x,y) be in IM. Thus 

1. If l{y) < l{x)^ there is an unique A G IM \ M such that X = (^,0), 



2. Ifl{y) > l{x), there is an unique AelR\R such that X = (0,A) = \(A, 0), 



3. If l{y) = l{x)^ there is an unique A = a eM. such that X = (a, 0) = (0, —a). 



Any element X — [A, 0) with A G IM — M is said positive and we write X > 0. Any element X — (0, A) with 
yl G M - K is said negative and we write X <0. We w rite X > A" if A" \ A" > 0. For example if X and X' 
are positive, X > X' <^==> 1{X) > 1{X').. The elements (a,0) with a G M* are neither positive nor negative. 

Ir^, one defines on the abelian group M , a structure of quasi linear space. Our approach is a little bit 
different. We propose to construct a real vector space structure. We consider the external multiplication: 

• :IRxill — >m. 

defined, for all A G IM, by 



a- (A,0) = {aA,0) 
a- {0,A) = {0,aA) 



for all a > 0. If a < we put /3 = —a. So we put: 



a-{A,0) ^ (0,/3A) , 
a-{0,A) -(/3A,0). 



We denote aX instead oi a ■ X. This operation satisfies 
1. For any a G M and A G IM we have: 



a{\X) = \{aX), 
{-a)X = \laX). 



2. For all a, /3 G M, and for all X, X' G M, we have 



(a + ;3)A' = aA" + ^A', 
a(A + A") =aA + aA", 
{a0)X ^a{/3X). 



Theorem 1 The triplet (IM, +,•) is a real vector space and the vectors X\ = ([0,1], 0) and Xi — ([1,1], 0) 
of IM determine a basis of IM. So dimR IM = 2. 

Proof. We have the following decompositions: 



([a,6],0) = (6 - a)Ai + 0X2, 
{0,[c,d]) = {c-d)Xi-cX2. 
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The linear map 

ifi-.m — >R^ 

defined by 

ip{ l[a,blO) ) = {b-a,a), 
^{{0,[c,d]))^{c-d,-c) 

is a linear isomorphism and IM is canonically isomorphic to M^. 

Remark. Let E be the subspace generated by A2. The vectors of E correspond to the elements which have 
a non defined sign. Then the relation < defined in the paragraph 1.2 gives an order relation on the quotient 
space m/E. 

C. A Banach structure on IR 



Any element <Y G IM is written {A, 0) or (0, A). We define its length 1{X) as the length of A and its center 
as c{A) or —c{A) in the second case. 

Theorem 2 The map \\ \\ : IM — > M. given by 

\\X\\=l{X) + \ciX)\ 
for any X Cz IR is a norm. 

Proof. We have to verify the following axioms: 

1) IJA-II =0-^=^ ^- = 0, 

2) VAeM \\XX\\ = \\\\\X\\, 

3) \\X + X'\\<\\X\\ + \\X'\\. 

1) If 11-^11 = 0, then 1{X) = \c{X)\ = and A" = 0. 

2) Let A e M. We have 

IIAA-II - 1{XX) + \c{XX)\ = \X\l{X) + \X\\c{X)\ - lAHlA-H. 



3) We consider that I refers to X and J refers to X' thus X = (1, 0) or = (0, /). We have to study the two 
different cases: 



i) If A" + A" = (/ + J, 0) or (0, 1 + J), then 

\\X + X'\\= 1(1 + J) + \c{I + J)\ = 1{I) + 1{J) + \c{I) + c{J)\ < 1{I) + \c{I)\ + 1{J) + |c(J)| 
= 11^11 + IIA'II. 



u) Let X + X' = (/, J). If (/, J) = {K, 0) then X + J = / and 

\\X + X'\\^ IIP^II ^ 1{K) + \c{K)\ = 1{I) 1{J) + |c(/) - c(J)| 
that is 

\\X + X'\\< 1(1) + |c(/)| - /(J) + |c(J)| < /(/) + |c(/)| + liJ) + |c(J)| = \\X\\ + \\X'\\. 
So we have a norm on IM. 
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Theorem 3 The normed vector space IM is a Banach space. 

Proof. In fact, all the norms on M^ are equivalent and M^ is a Banach space for any norm. The vector space 
IM is isomorphic to M^. Thus it is complete. 

Remarks. 

1. To define the topology of the normed space IM, it is sufficient to describe the e-neighborhood of 
any point xo G IK for e a positive infinitesimal number. We can give a geometrical representation, 
considering xo = ([a, &],0) represented by the point (a, 6) G M? . We assume that xo — ([a, ^],0) and e 
an infinitesimal real number. Let ^i, • • • , ^4 the points Ai — {a — e,h — e), A2 = (a + %,b — |), ^3 — 
(a + £,& + e), A4 = (a— §,&+ |). If < a < h, then the e-neighborhood of xo — ([a, 6],0) is represented 
by the parallelograms whose vertices arc yli, A2, A3, ^4. 

2. We can consider another equivalent norms on IM. For example 

\\X\\^\\-^X\\^Sup{\xl\y\) 
where X = ([a;,y],0). But we prefer the initial one because it has a better geometrical interpretation. 

III. DIFFERENTIAL CALCULUS ON II 



As M is a Banach space, we can describe a notion of differential function on it. Consider Xq = (Xq, 0) 
in IM . The norm ||.|| defines a topology on IM whose a basis of neighborhoods is given by the balls 
B{XQ,e) = {X e M, IIA" \ A'oll < e}- Let us characterize the elements oi B(XQ,e). X^ = (Xo,0) = ([a, 6],0). 



Proposition 2 Consider Xq = {Xq,0) = ([a, fe],0) in IM and e ~ 0, £ > 0. Then every element of B{Xo,s) 
is of type X = {X, 0) and satisfies 

1{X) e Br(?(Xo),£i) andc{X) G Sr(c(Xo),£2) 

with £1, £2 > and £1 + £2 < £, where B^{x, a) is the canonical open hall in M of center x and radius a. 

Proof. First case : Assume that X = (AT, 0) — {[x,y],0) . We have 



X\Xo = {X,Xo)^i[x,y],[a,b]) 

^ r {[x-a,y-b],0) if 1{X) > liXo) 
1 {0,[a-x,b~y])iil{X)<l{Xo) 

If 1{X) >liXa) we have 

y — b + X — a 



\X\Xo\\ = iy-b)-{x-a) + 



2 

= 1{X) ~ liXo) + \c{X) - c{Xo)\. 

As 1{X) - 1{Xq) > and \c{X) - c{Xq)\ > 0, each one of this term if less than s. If 1{X) < /(Aq) we have 

\\X \ A-oll = l{Xo) - 1{X) + |c(Ao) - c(A)|. 
and we have the same result. 



Second case : Consider X = (0, A) = ([x,j/],0) . We have 



X\Xo = (0, Ao + A) = {[x + a,y + b]) 
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and 

\\X \ XoW - l{Xo) + 1{X) + |c(Xo) + c{X)\. 
In this case, we cannot have \\X \ Ao|| < e thus X ^ ,B(Xo,£). 

Definition 4 j4 Junction f : M — ^ M is continuous at A'o «/ 

Ve > 0, 377 > sMc/i t/iaf HA" \ A'o|| < 77 implies \\.f{X) \ /(A'o)ll < e. 

Consider (^"1, ^"2) the basis of IM given in section 2. We have 

f{X) = /i (<¥)<¥! + /2(A')A'2 with fr.m^R. 

If / is continuous at Xq so 

/(A-) \ f{Xo) - (/il-Y) - .hiX„))Xi + ihiX) - f2{Xo))X2. 

To simplify notations let a = fi{X) - fi{Xo) and /3 ^f2{X) - f2{Xo). If \\f{X) \ /(-Yo)|| < e, and if we 
assume fi{X) — /i(A'o) > and j^iX) — /2(<^o) > (other cases are similar), then we have 



;(aA'i+;3A'2) = ?([/3,a + /3],0) <e 
thus fx{X) - /i(A:'o) < e. Similarly, 



c(aA'i+/3A'2)==c([/3,a + /3],0)= 3 + /^ < ^ 
and this implies that /2('^) ^ /2('^o) < £• 
Corollary 5 f is continuous at Xq if and only if fi and f^ are continuous at Xq. 

Definition 6 Consider Xq in IM and f : IM — > IM continuous. We say that f is differentiable at Xq if 
there is g : IM — ;■ IM linear such as 

\\f{X) \ f{Xo) \ g{X \ Xo)\\ - o{\\X \ A^oll). 
Examples. 

• f{X) — X. This function is continuous at any point and differentiable. It's derivative is /'(<-f) = 1. 



• f{X) = X^. Consider Xq = (Xo,0) = ([a,6],0) and X e B{Xo,e). We have 

\\X^ \ X^W ^ \\{X \ Xq){X + Xo)\\ 
< \\X \ Xq\\\\X + Xo\\. 

Given e > 0, let ?7 = t—- -— r, thus if jjA:' \ Ao|| < 77, we have \\X'^ \ Xq\\ < e and / is continuous 

and differentiable. It is easy to prove that f'{X) = 2X is its derivative. 

• Consider P ^ qq + aiX + h a„X" e M[X]. We define / : !M — ^ IM with f{X) =- 00-^2 + aiX + 

• • • + a^X" where X" = X ■ X"^^ . From the previous example, all monomials are continuous and 
differentiable, it implies that / is continuous and differentiable as well. 

• Consider the function Q2 given by Q2{[x,y]) — [x'^,y'^] if |x| < \y\ and Q2{[x,y] — [y'^,x'^] in the other 
case. This function is not differentiable. 
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IV. A 4-DIMENSIONAL ASSOCIATIVE ALGEBRA ASSOCIATED WITH II 

In introduction, we have observed that the semi-group IM is identified to Pi.i U7'i,2 UPi.a. Let us consider 
the following vectors of M^ 

ei-(l,l), 
62 -(0,1), 

63 = (-1,0), 

64 = (-1,-1). 

They correspond to the intervals [1, 1], [0, 1], [—1, 0], [—1, —1]. Any point of T'l^i U 7^1,2 U V1.3 admits the 
decomposition 



(a, b) — aiei + 0:262 + 0:363 + 0:464 
with Oi > O.The dependance relations between the vectors Ci are 

62 = 63 + 61 

64 = -61. 

Thus there exists a unique decomposition of (a, b) in a chosen basis such that the coefhcicnts are non negative. 
These basis are {61,62} for T'l,!, {62, 63} for 7^1,2, {63, 64} for 7^1,3, Let us consider the free algebra of basis 
{61,62,63,64} whose products correspond to the Minkowski products. The multiplication table is 





61 


62 


63 


64 


ei 


61 


62 


63 


64 


62 


62 


62 


63 


63 


63 


63 


63 


62 


62 


64 


64 


63 


62 


61 



This algebra is associative. Let (p 
its image in Ai, we have: 



— > Ai the natural injective embedding. If we identify an interval with 
Theorem 7 The multiplication of intervals in the algebra A4 is distributive with respect the addition. 



The application is not bijective. Its image on the elements X — (x, 0) = ([a, 6], 0) is: 

X = [a, b] e Vi,i,(p{X) = a6i + (6 - a)e2 (a > 0, 6 - a > 0) 

X = [a, b] e ri^2,(pW) = -a63 + 662 (-a > 0, 5 > 0) 

X = [a, b] e ^^1^3, ifilx) = -664 + (6 - 0)63 (-fe > 0, 6 - a > 0). 

Consider in A4 the linear subspace F generated by the vectors 61 — 62 + 63, 61 + 64. As 

(61 + 64)(6i + 64) = 2(61 + 64) 
(61 + 64) (61 - 62 + 63) = 61 + 64 
(61 - 62 + 63)(6i - 62 + 63) = 61, 

F is not a subalgebra of A^. Let us consider the map 

^ : m -^Ai/F 

defined from (p and the canonical projection on the quotient vector space A^/ F. A vector x = ^ oiiCi G .44 
is equivalent to a vector of Aa with positive components if and only if 

02 + 03 > 0. 
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In this case, all the vectors equivalent to a; = ^ a^ei with a2 + a^ > correspond to the interval [ai — 
a^ — a4, ai + a2 ~ 04] of IM. Thus we have for any equivalent classes of A4/F associated with ^ UiCi with 
Q^2 + c^a ^ a prcimage in IM. The map ^ is injective. In fact, two intervals belonging to pieces 'Pi.i,'Pij- 
with i y^ j, have distinguish images. Now if (a, b) and (c, d) belong to the same piece, for example Pi,i, thus 

(/3(a, b) = {(a + A + /i, 6 — a — A, A, /i). A, /i G M.} 

If Tp(c, d) = ^(a, b), there are A, /i G M such that {c,d) = (a + X + fi,b — a — X, X, fi). This gives a — c,b ~ d. 
We have the same results for all the other pieces. Thus Tp : IM -^Aa/F is bijective on its image, that is the 
hyperplane of A4/F corresponding to a2 + 0:3 > 0. 

Practically the multiplication of two intervals will so be made: let X,Y eR. Thus X — J2 oa^i, ^ = XI A^i 
with ai,/3j > and wc have the product 



this product is well defined because (p{X).ip{Y) G Inup. This product is distributive because 

X*{Y + Z) =i^-\>f{X).if{Y + Z)) 

= W\^{X).p{Y)+^{X).^{Z)) 
= X»Y + X»Z 



Remark. We have 



W\v{X).v{Y + Z)) ^ ^-'{^{X)).^~\p{Y + Z))). 



We shall be careful not to return in IM during the calculations as long as the result is not found. Otherwise 
wc find the semantic problems of the distributivity. 

We extend naturally the map (p : IM -^Aa to IM by 



^(AO) ^ ^{A) 
^(0,A) = -^(A) 

for every ^ G IM. 

Theorem 8 The multiplication 

x'mx" = ip-\ip{x').ip{A:")) 

is distributive with respect the addition. 

Proof. This is a direct consequence of the previous computations. 



In A4 we consider the change of basis 



ei -62 



6^ — 6^ , Z — ^, ^ 
64 = 64 - 63. 

This change of basis shows that A4 is isomorphic to A4 





61 


62 


63 


64 


ei 


61 








64 


62 





62 


63 





63 





63 


62 





64 


64 








61 



The unit of A'4 is the vector ei +62. This algebra is a direct sum of two ideals: A'4 ^ h + I2 where /i 
is generated by 61 and 64 and I2 is generated by 62 and 63. It is not an integral domain, that is, we have 
divisors of 0. For example 61 • 62 = 0. 
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Proposition 3 The multiplicative group A^ of invertible elements of A4 is the set of elements x — 
{xi,X2,X3,X4) such that 

X4 ^ ±Xi, 
^ X3 ^ ±X2- 

If X ^ AX we have: 

-if Xi X2 X3 X4^ 



X 



Jl "^ A t<J Q iXj O i/v Q i/v Q ^1 -^ A 



Let us compute the product of intervals using the product in Ai and we compare with the Minkowski 
product. Let X = [a, b] and Y = [c, d] two intervals. 

Lemma 1 If X and Y are not in the same piece Vi^i, then X uY corresponds to the Minkowski product. 

Proof, i) If X e Vis and Y e Vi^ then ip{X) = {a,b- a, 0, 0) and (f{Y) = (0, d, -c, 0). Thus 

ip{X)ip{Y) = (aei + (6- a)e2)(de2 - ces) 
= bde2 — cbes 
= (0, bd, -cb, 0) 
= ^{[cb,bd]). 

ii) If X e Vi,! and Y e Vi^s then Lp{X) = {a,b- a, 0, 0) and (f{Y) ^ (0, 0,d-c, -d). Thus 

ip{X)ip{Y) = {aei + {b- a)e2){{d - c)e3 - de^) 
= [ad — bc)e-^ — ade^ 
— (0, 0, ad — cb, —ad) 
= (p{[bc,ad]). 

iii) If X e Pi,2 and Y e Pi,3 then (p{X) = (0, b, -a, 0) and Lp{Y) = (0, 0, rf - c, -d). Thus 

ip{X)ip{Y) = {be2 - ae3){{d - c)e3 - dei) 
= ace2 — bce^ 
= (0, ac, -cb, 0) 
= (p{[bc,ad]). 

Lemma 2 If X an Y are both in the same piece T'l.i or Vi^3, then the product X uY corresponds to the 
Minkowski product. 

The proof is analogous to the previous. 

Let us assume that X — [a,b] and Y — [c,d] belong to Pi. 2- Thus (p{X) — (0, 6, — a, 0) and ip{Y) = 
(0,d,-c,0). We obtain 

XY = (6e2 — 063) ((ie2 — ce^) = {bd + ac)e2 + {—be — ad)e3. 

Thus 

[a, b] [c, d] = [be + ad, bd + ac] . 

This result is greater that all the possible results associated with the Minkowski product. However, we have 
the following property: 

Proposition 4 Monotony property: Let Afi, ^"2 G IK. Then 

XiCX2=^ Xi» Z CX2» Z for all Z G M. 

^{Xi) < ^{X2) =^ ^{Xi . Z) < ^{X2 . Z) 
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The order relation on Ai that ones uses here is 



r (xi,a;2,0,0)<(yi,2/2,0,0) 
(a^i, 2^2, 0,0)<(0,y2, 2/3,0) 
(0,a;2,a:3,0)< (0,^2,2/3,0) 
(0,0,a;3,a:4) < (0,^2,2/3,0) 

I (0,0,a;3,X4)< (0,0,2/3,2/4) 



yi < xi and X2 < 2/2, 

X2 < 2/2, 

X3 < 2/3 and X2 < 2/2, 
2^3 < 2/3, 
X3 < 2/3 and y4 < 2:4. 



Proof. Let us note that the second property is equivalent to the first. It is its translation in Ai- We can 
suppose that Xi and X2 sue intervals belonging moreover to 7^1.2: <p(<^i) = (0, b, —a, 0), (f{X2) = (0, d, — c, 0). 

li (p{Z) = (21,22,2:3,^4), then 

^{Xi • Z) — (0, bzi + bz2 — az3 — 024, —azi + bz^ — az2 + 624, 0), 
(^(^2 ' Z) ~ (0, dzi + dz2 — CZ3 — cz/^, —czi + dz3 — CZ2 + dz^, 0). 



Thus 



(fe - d){zi + Z2) - (a - c)(z3 - Z4) < 0, 
-(a - c)(zi + Z2) + {b- d){z3 = Z4) < 0. 



^{X, • Z) < ip{X2 . Z) ^ 
But (6 — d), — (a — c) < and 22,23 > 0. This implies ^{Xi • Z) < Tp(X2 • Z 



V. THE ALGEBRAS A„ AND AN BETTER RESULT OF THE PRODUCT 



We can refine our result of the product to come closer to the result of Minkowski. Consider the one 
dimensional extension ^40lRe5 = ^5, where 65 is a vector corresponding to the interval [—1, 1] of T^i 2. The 
multiplication table of A5 is 





61 


62 


63 


64 


65 


ei 


61 


62 


63 


64 


65 


62 


62 


62 


63 


63 


65 


63 


63 


63 


62 


62 


65 


64 


64 


63 


62 


61 


65 


65 


65 


65 


65 


65 


65 



The piece Pi_2 is written 7^1,2 = 'Pi, 2,1 U Pi, 2,1 where 7^1,2.1 = {[a, b],—a < b} and 7^1,2,2 = {[a, b], —a > b}. 
liX = [a, b] € 7^1,2,1 and Y = [c, d\ G Vi,2,2,'t'hus 

ip{X).Lp{Y) = (0, fe + a, 0, 0, -a).(0, 0, -c - d, 0, d) = (0, -(a + 6)(c + d), 0, 0, a(c + d) + bd). 

Thus we have 

X uY = [—bd — ac — ad, —6c]. 



Example Let X = [-2, 3] and Y ^ [- 



The product in A^ gives 



The Minkowski product is 



-4, 2]. We have X e 7^1,2,1 and Y e 7^1,2,2 • The product in A4 gives 
X.r = [-16,14]. 

X»Y = [-12,10]. 



[-2,3].[-4,2] = [-12,i 
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Thus the product in A5 is better. 

Conclusion. Considering a partition of 7^i,2, we can define an extension of ^4 of dimension n, the choice 
of n depends on the approach wanted of the Minkowski product. For example, let us consider the vector eg 
corresponding to the interval [—1, 2]- Thus the Minkowsky product gives eg. eg = 67 where 67 corresponds 
to [— i , 1] . We obtain a 7-dimensional associative algebra whose table of multiplication is 





61 


62 


63 


64 


65 


66 


67 


ei 


61 


62 


63 


64 


65 


66 


67 


62 


62 


62 


63 


63 


65 


66 


67 


63 


63 


63 


62 


62 


65 


67 


66 


64 


64 


63 


62 


61 


65 


67 


66 


65 


65 


65 


65 


65 


65 


65 


65 


66 


66 


66 


67 


67 


65 


67 


66 


67 


67 


67 


66 


66 


65 


66 


67 



Example Let X 

coefhcients writes 



-2,3] and Y = [—4,2]. The decomposition on the basis {ei,--- ,67} with positive 



Thus 



X = 65 + 267, r = 266 



X.r=(65 + 267)(466)=465 



'12,81 



We obtain now the Minkowski product. In general, when one increases the algebra dimension, the product 
will be closer to the Minkowski one and one still get the distributivity and associativity. 



VI. NUMERICAL IMPLEMENTATION 

In this section, we show some examples of interval arithmetics applications on simple problems which will 
prove how this new approach efficient and robust is. 



A. Arithmetic implementation in python 

We have choosen python programming langage^. The main reason is that it is a free object-oriented 
langage, with a huge number of numerical libraries. One of the main advantage of python is that first it is 
possible to link the source code with others written in C/C++, FORTRAN, and second, it interacts easily 
with other calculations tools such as SAGE?^ and Maxima^*^ in order to do formal calculations with python 
langage. But here, we present pure numerical applications within python environnement. The translation in 
other langages such as C++ and scilao^ is very easy and would be available soon. To start it is necessary 
to import the intervalJib library which has been developped to define intervals, vector and matrices of 
intervals, and all the arithmetic operations. 

Python 2.6.6 (r266: 84292, Sep 15 2010, 15:52:39) 

[GCC 4.4.5] on linux2 

Type "help", "copyright", "credits" or "license" for more information. 

>» from interval_lib import * 



The instanciation of an interval [x,y] is done with interval{x,y, order). The variable order corresponds to 
the dimension of the algebra used to represent the intervals. Its value is set to 4 by default, which is the 
minimal one. Another way to define an interval such as [x — e,x + e] is inter val{x,eps — e). Now, let's define 
the intervals [—1,2], [3,4], [3, 12] and [1,3] for example : 
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>» a=interval(-l,2) 
>» b=interval(3,4) 
>» c=interval(3,12) 
>» d=interval(2,eps=l) 

It is possible to have more information on each interval : 

>» print c.min,c.max 

3.0 12.0 

>» print abs(c) , c .width.c .midpoint 

16.5 9.0 7.5 

A partial order relation can be implemented on the set of intervals IM : 

Vx, y £ IM, X (t y, center{x) < center{y) ^ x < y (1) 

and 

Vx, y G IM, X <Z y, width{x) < width{y) <^ x < y (2) 

This can be extended to a total order relation. 



>» print a,b,a<b 
[-1.0,2.0] [3.0,4.0] True 

>» print c,d,d<c 
[3.0,12.0] [1.0,3.0] True 



B. Semantic and True Arithmetic 

There are two possible arithmetic implementations, depending on the choosen semanticP^Hl, in the first 
one, called semantic arithmetics, the substraction of two intervals x and y is done according io x — y — x -\- 
(-y), and the addition of those two terms. For example [2,3]-[0,l] = [2,3] + [-l,0] = [2+(-l),3+0] = [1,3], 
and [— 1, 1] — [— 1, 1] — [— 1, l] + [— 1, 1] = [—2,2] ^ [0,0] — 0. As mentionned in the introduction of this paper, 
"negative" intervals do not have a physical meaning and the addition/substraction between two intervals 
can not be easily transfered to the bounds of the resulting interval. This yields to the fact that differential 
calculus in this framework is not relevant and one has to compute the derivatives in the center of the intervals 
in order to recover a certain meaning. It is not obvious to transfer natural functions to inclusion ones. In 
the second framework, called true arithmetic, the substractions are done in the algebra An with n > A. 
For the previous example in An : [2,3] - [0,1] = (2,1,0,0) - (0,1,0,0) = (2,0,0,0) = [2,2] and for any 
interval x, x — x = [Q,Q\ = 0. In this arithmetic, it is possible to perform differential calculus and to transfer 
natural functions to inclusion ones by replacing the terms in the definition by intervals. One has to note 
that in both cases, multiplication remains distributive and associative according to addition. But division is 
distributive for substraction only for the true arithmetic even if it is distributive for addition in the semantic 
one. The main reason of this phenomenon, is that the opposite intervals have no real meaning, and it 
remains to the user to modelize correctly the physical problem. Moreover, there is no wrapping effects and 
data dependencies as shown on simple examples below. 
Some other examples : 

>» # Semantic arithmetic 

>» print a-a,a*b,b*a,b/b,c+l 

[-3.0,3.0] [-4.0,8.0] [-4.0,8.0] [1.0,1.0] [4.0,13.0] 

and 
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>» # True arithmetic 

>» print a-a,a*b,b*a,b/b,c+l 

[0.0,0.0] [-4.0,8.0] [-4.0,8.0] [1.0,1.0] [4.0,13.0] 

The division is not allowed for intervals containing : 

>» print b/a 

Interval division in A4 not allowed ! ! 

Let's see the distributive operations : 

>» # True arithmetic 

>» print a* (b+c) , a*b+a*c , (a+b) /c , a/c+b/c 

[-16.0,32.0] [-16.0,32.0] [0.5,0.916666666667] [0.5,0.916666666667] 

>» print a*(b-c) ,a*b-a*c, (a-b)/c,a/c-b/c 

[-16.0,8.0] [-16.0,8.0] [-1.08333333333,-0.166666666667] [-1.08333333333,-0.166666666667] 

In the semantic arithmetic 

# Semantic arithmetic 

>» print a* (b+c) , a*b+a*c , (a+b) /c , a/c+b/c 

[-16.0,32.0] [-16.0,32.0] [0.5,0.916666666667] [0.5,0.916666666667] 

>» print a* (b-c) , a*b-a*c , (a-b) /c , a/c-b/c 

[-28.0,20.0] [-28.0,20.0] [-0.833333333333,-0.416666666667] [-1 .08333333333,-0 . 166666666667] | 

In the semantic framework, the distributivity of division according substraction is lost but not according 

addition. This is due to the calculation of a — 6 before to be divided by c. However the division distributivity 

is always fully respected in the true arithmetic. 

Another interesting example shows that one gets no wrapping and data dependancy for the two arithmetic 

frameworks. 

>» def fl(x):return x**2-2*x+l 
>» def f 2 (x): return x*(x-2)+l 
>» def f3(x):return (x-l)**2 

>» # Semantic arithmetic 
>» print fl(a), f2(a), f3(a) 
[-7.0,8.0] [-7.0,8.0] [-7.0,8.0] 
>» print fl(b), f2(b), f3(b) 
[2.0,11.0] [2.0,11.0] [2.0,11.0] 

and 

>» # True arithmetic 
>» print fl(a), f2(a), f3(a) 
[-1.0,2.0] [-1.0,2.0] [-1.0,2.0] 
>» print fl(b), f2(b), f3(b) 
[4.0,9.0] [4.0,9.0] [4.0,9.0] 

One remarks that the true arithmetic results are always included in the ones obtained with the semantic 
arithmetic. 
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C. Optimization examples 



1. Fixed-step gradient descent method 



Here is a script example of minimization with fixed-step gradient method which belongs to the so-called 
gradient descent method^. This algorithm and this example are very simple but it shows that the result is 
garanted to be found within the final interval. 



Fixed step gradient optimization of x*exp(x| 




Lower bound - 

Center of interval - 

Upper bound - 



500 1000 1500 2000 2500 3000 3500 
Iterations 



Figure 1. Convergence of the fixed-step gradient algorithm for the function x i-^ x ■ exp(x) to an interval centered 
around —1. 



from interval_lib import * 

# Example of fixed step gradient descent method 

f ile=open("res .data" , "w") # Data file to be plotted 

h=l.e-6 # Finite difference step 

def f(x):return x*(exp(x)) # Function to be minimized 

def fp(x):return interval( ( (f (x+h)-f (x-h) ) .midpoint)/h/2. ) # Finite difference 

x=interval(2, eps=.l) # Initial guess 

rho=interval(l . e-2) # Gradient step 

epsilon=l . e-6 # Accuracy of the gradient 

while abs(f p(x) )>epsilon: # Descent loop 

fprime=f p(x) 

x=x-rho*f prime 

f ile. write (("7,f 'Li °Li °Lf Zf \n")7o(x.min, x.max, x. midpoint, fprime.min, fprime.max)) 
file . closeO 

In the true arithmetic, the finite differences are "smaller" and it has meaning to do derivative calculations. 
This is due to the fact that for close intervals, the difference is close to 0. One has just to change 

def fp(x):return (f (x+h)-f (x-h) )/h/2. # Finite difference 



The result shown in figure [2] is impressive, because for any initial guess the interval width decreases to 
converge to real point minimum. In the semantic interval on figure [11 the width of the interval does not 
decrease and the center converges to the right value. This is due to finite difference calculation at the center. 
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Optimization of x*exp(x) with fixed-step gradient method 




lOOOC 



Figure 2. Convergence of the fixed-step gradient algorithm with true arithmetic for the function x t~^ x ■ exp(a::) to 
an interval centered around —1. 



2. Newton-Raphson method 



Let's optimize the same function x t-^ x ■ exp(a;) with a second order method such as the Newton-Raphson 
one, which is the basis of all second order methods such as Newton or quasi-Newton's oneg^^. It finds the 
same minimum which is an interval centered around —1. 



# Example of Newton-Raphson method 

from interval_lib import * 

f ile=open("res .data" , "w") # Data file to be plotted 

h.= (l.e-6) # Finite difference step 

def f(x):return x*(exp(x)) # Function to be minimized 

def fp(x):return intervalC ( (f (x+h)-f (x-h) ) .midpoint)/h/2. ) # Finite difference 

def fp2(x) :return intervalC ( (f(x+h)+f (x-h) -2*f (x) ) .midpoint) /(li*h) ) # Finite difference 

x=interval(2, eps=.l) # Initial guess 

epsilon=l . e-10 # Accuracy of the gradient 

while abs(f p(x) )>epsilon: # Descent loop 

fprime=f p(x) 

f second=f p2 (x) 

file .write ( ("7of "/of 7„f 7,f 7of \n")yo(x.min, x.max, x. midpoint, fprime.min, fprime.max)) 

x=x-fprime/f second 
file . closeO 



Another interesting example is shown on the figures 5[6 and [7] for different initial guess intervals. We 
would like to optimize x i— >■ (a;^ — 1)^- One has to change the finite differences calculated in the center of the 
interval by classical finite differences : 



def fp(x) : return (f (x+h)-f (x-h) )/h/2. # Finite difference 

def fp2(x) : return (f (x+h)+f (x-h)-2*f (x) )/(h*h) # Finite difference 



The minima are {—1, 0, 1}. One can see that depending on the initial guess, this simple algorithm finds the 
right real point minima. 
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Newton-Raphson optimization of x*exp(x) 




Lower bound - 
Center of interval - 
Upper bound - 



12 3 4 5 6 7 8 

Iterations 

Figure 3. Convergence of the Newton-Raphson algorithm for the function x t-^ x ■ exp(a;) to an interval centered 
around —1. 



Optimization of x*exp(x) with Newton-Raphson method 




Iterations 



Figure 4. Convergence of the Newton-Raphson algorithm with true arithmetic for the function x i~^ x ■ exp(a::) to an 
interval centered around —1. 



D. Matrix diagonalization and inversion 



1. Diagonalization 



As an example, we define the matrix M whose elements are intervals centered around a certain real number 
with a radius e. 



M 



[l-e,l- 
[3-6,3 



e] [2-e,2- 
el [4-e,4- 



If one uses scilah to compute the spectrum of the previous matrix without radius (e = 0), the highest 
eigenvalue is approximatively 5.3722813 and the corresponding eigenvector is (0.4159736,0.9093767). In 
order to show that arithmetics and interval algebra developped above is robust and stable, let's try to 
compute the highest eigenvalue of an interval matrix. One uses here the iterate power method, which is very 
simple and constitute the basis of several powerful methods such as deflation and others. The two figures 
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Minimization of (x**2-l)**2 witli Newton-Raphson mettiod 




200 400 



1000 1200 1400 1600 



Iterations 



Figure 5. Convergence of the Newton-Raphson algorithm with true arithmetic for the function x i 



Minimization of (x**2-l)**2 witli Newton-Raphson method 



-0,5 r 



^ -1,5 



c -2,5 



•3 - 



-3,5 L 



12 3 4 5 6 7 

Iterations 



Figure 6. Convergence of the Newton-Raphson algorithm with true arithmetic for the function x >->■ {x^ — 1)^. 

[8] and [9] show clearly for different value of e the stability of the multiplication, and the largest eigenmode 
is recovered when e = 0. The other eigen modes can be computed with the deflation methods for example 
which consists to withdraw the direction spanned by the eigenvector associated to the highest eigenvalue to 
the matrix by constructing its projector and to do the same. Several methods are available and efficient. 
We have choosen to compute only the highest eigenvalue and its corresponding eigenvector in order to show 
simply the efficiency of our new artihmetic. The corresponding code in python is described below : 



# Example of an interval matrix diagonalization 
from interval_lib import * 

f ile=open("res .data" , "w") # Data file to be plotted 
for i in xrajige(lO) :# Loop on the radius of the matrix elements 
epsilon=10 . ** (-i) 

# Construction of the matrix 

a=interval(l , eps=epsilon) ;b=interval(2, eps=epsilon) 
c=interval(3, eps=epsilon) ;d=interval(4, eps=epsilon) 
u=Vector([a, b] ) ; v=Vector ( [c, d] ) 
uO=Vector( [interval (1) , interval(l)] ) #lnitial guess 
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Minimization of (x**2-l)**2 with Newton-Raphson method 



3,5 r 
3 
^ 2,5 

C 
J 



^ 2 
> 

£ 1.5 : 

1 ■ 

0,5 



12 3 4 5 6 7 

Iterations 



Figure 7. Convergence of the Newton-Raphson algorithm with true arithmetic for the function x i 



Iterate power method 



"^ 3,5 
■f le- 




le-05 0,0001 0,001 0,01 
intervai radius of matrix eiements 



Figure 8. Comparison of the largest eigenvalue found with scilab 5.3722813 versus the iterate power method one 
computed with intervals. 



m=Matrix( [u, v] ) # Interval matrix to be diagonalized 
e,v=iterate_power (m, uO, 10) #Power iteration 

file .write ("7of Zf °/f 7of "/of 7of %f \n"7o(epsilon, e.min, e.max, v[0].min, v[0] .max,v[l] .min, v[l].max ))| 
file . closeO 



2. Inversion 



Let's define a symmetric matrix. 




[5 - e, 5 + e] 
[6 - e, 6 + e] 
[3 - e, 3 + e] 
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Iterate power method 




0,2 
le- 



le-05 0,0001 0,001 0.01 
Interval radius of matrix elements 



Figure 9. Comparison of the eigenvector corresponding to the largest eigenvalue found with scilab 

(0.4159736, 0.9093767) versus the iterate power method one computed with intervals. 



We would like to use the well-know Schutz-Hotelling algorithnF^*^ to inverse a matrix X 



Xo = ^J^, X, = X,_i(2 - A ■ X,_i), Vn > 1 



(3) 



Scilab gives numerically for e = 



-0.2678571 0.1607143 0.125 
M-^ = I 0.1607143 -0.1964286 0.125 



0.125 



0.125 -0.125 



(4) 



The python code is very simple : 



# Example of an interval matrix inversion 
from interval_lib import * 

epsilon=.2 #intervals radius 

a=interval(l, eps=epsilon) ;b=interval(2, eps=epsilon) ; c=interval(3, eps=epsilon) 

d=interval(4, eps=epsilon) ;e=interval(5, eps=epsilon) ;f =interval(6, eps=epsilon) 

# Build the matrix m 

u=Vector([a, d, e] ) ; v=Vector ( [d, b, f ] ) ;w=Vector ( [e , f, c] ) 
m=Matrix([u, v, w] ) ; inv_m=schultz(m) ; inv_inv_m=schultz(inv_m) 

# Display results 
print "M=", m 

print "Inverse matrix = ", inv_m 

print "M'(-1)*M=", inv_m*m 

print "M*M~(-1)=", m*inv_m 

print "(M~(-l) )"(-!)=" ,inv_inv_m 

We obtain with intervals for e = 0.2 for example : 

M= [* 

[0.8,1-2] [3.8,4.2] [4.8,5.2] 

[3.8,4.2] [1.8,2.2] [5.8,6.2] 

[4.8,5.2] [5.8,6.2] [2.8,3.2] 

*] 
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Inverse matrix = [* 

[-0 . 267918088737 , -0 . 267790262172] [0 . 160409556314 , . 161048689139] [0 . 12457337884 , . 125468164794] | 

[0 . 160409556314 , . 161048689139] [-0 . 19795221843 , -0 . 194756554307] [0 . 122866894198 , . 12734082397] | 

[0 . 12457337884 , . 125468164794] [0 . 122866894198 , . 12734082397] [-0 . 127986348123 , -0 . 121722846442] | 

*] 

M~(-1)*M= [* 

[1 . , 1 • 0] [-2 . 77555756156e-17 , . 0] [-1 . 80411241502e-16 , -1 . 66533453694e-16] 

[-2 . 22044604925e-16 , -2 . 11636264069e-16] [1 . , 1 • 0] [-1 . 38777878078e-17 , . 0] 

[0.0,1- 38777878078e-17] [0 . , . 0] [1 . , 1 . 0] 

*] 

M*M'(-1)= [* 

[1 . , 1 . 0] [-1 . 21430643318e-16 , -1 . 11022302463e-16] [0.0,1. 38777878078e-17] 

[0.0,0.0] [1.0,1.0] [-1.38777878078e-17,0.0] 

[4 . 16333634234e-17 , 5 . 55111512313e-17] [-1 . 80411241502e-16 , -1 . 66533453694e-16] [1.0,1.0] 

*] 

(M-(-l))-(-l)= [* 

[0.8,1.2] [3.8,4.2] [4.8,5.2] 

[3.8,4.2] [1.8,2.2] [5.8,6.2] 

[4.8,5.2] [5.8,6.2] [2.8,3.2] 

*] 

and for e = 0.1 

M= [* 

[0.9,1.1] [3.9,4.1] [4.9,5.1] 

[3.9,4.1] [1.9,2.1] [5.9,6.1] 

[4.9,5.1] [5.9,6.1] [2.9,3.1] 

*] 

Inverse matrix = [* 

[-0 . 267888307155 , -0 . 267824497258] [0 . 160558464223 , . 160877513711] [0 . 124781849913 , . 125228519196] | 

[0 . 160558464223 , . 160877513711] [-0 . 197207678883 , -0 . 195612431444] [0 . 123909249564 , . 126142595978] | 

[0 . 124781849913 , . 125228519196] [0 . 123909249564 , . 126142595978] [-0 . 12652705061 1 , -0 . 123400365631] | 

*] 

M~(-1)*M= [* 

[1 . , 1 . 0] [2 . 22044604925e-16 , 2 . 35922392733e-16] [1 . 66533453694e-16 , 1 . 7694179455e-16] 

[0.0,0.0] [1.0,1.0] [0.0,0.0] 

[-1 . 31838984174e-16 , -1 . 11022302463e-16] [-1 . 2490009027e-16 , -1 . 11022302463e-16] [1.0,1.0] 

*] 

M*M'(-1)= [* 

[1.0,1.0] [0.0,0.0] [-6.93889390391e-18,0.0] 

[0.0,0.0] [1.0,1.0] [0.0,6.93889390391e-18] 

[-3.46944695195e-18,0.0] [0.0,0.0] [1 .0, 1 .0] 

*] 

(M-(-l))-(-l)= [* 

[0.9,1.1] [3.9,4.1] [4.9,5.1] 

[3.9,4.1] [1.9,2.1] [5.9,6.1] 

[4.9,5.1] [5.9,6.1] [2.9,3.1] 

*] 

and for e — 0.01 

M= [* 

[0.99,1.01] [3.99,4.01] [4.99,5.01] 

[3.99,4.01] [1.99,2.01] [5.99,6.01] 

[4.99,5.01] [5.99,6.01] [2.99,3.01] 

*] 

Inverse matrix = [* 
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[-0 . 267860324247 , -0 . 267853946662] [0 . 160698378764 , . 160730266691] [0 . 124977730269 , . 125022373367] | 

[0 . 160698378764 , . 160730266691] [-0 . 196508106182 , -0 . 196348666547] [0 . 124888651345 , . 1251 1 1866834] | 

[0. 124977730269,0. 125022373367] [0 . 124888651345,0 . 125111866834] [-0. 125155888117,-0. 124843386433]| 

*] 

M~(-1)*M= [* 

[1 . , 1 . 0] [1 . 11022302463e-16 , 1 . 12323345069e-16] [5 . 55111512313e-17 , 5 . 63785129692e-17] 

[-1.517883041486-18,0.0] [1.0,1.0] [-5. 74627151417e-17,-5. 55111512313e-17] 

[-8 . 673617379886-19 , . 0] [-2 . 229119666636-16 , -2 . 22044604925e-16] [1.0,1.0] 

*] 

M*M"(-1)= [* 

[1 . , 1 . 0] [-2 . 168404344976-19 , . 0] [0 . , . 0] 

[-2.220446049256-16,-2.216109240566-16] [1.0,1.0] [0.0,0.0] 

[-1 . 123233450696-16 , -1 . 110223024636-16] [-5 . 57279916658e-17 , -5 . 551115123136-17] [1.0,1.0] 

*] 

(M-(-l))-(-l)= [* 

[0.99,1.01] [3.99,4.01] [4.99,5.01] 

[3.99,4.01] [1.99,2.01] [5.99,6.01] 

[4.99,5.01] [5.99,6.01] [2.99,3.01] 

*] 

It is obvious that this method is very stable and confirms that the true arithmetic operations are robust. It 
is not difficult to extend usual linear iterative algebra numerical algorithms to intervals. It permits to solve a 
lot problems where the entries of the matrices are not well defined, especially for automation applications^. 



VII. CONCLUSION 

We have presented a better algebraic way to do calculations on intervals. This approach"'^ is done by 
embedding the space of intervals into a free algebra of dimension greater or equal to 4. This permits to 
obtain all the basic arithmetic operators with distributivity and associativity. We have shown that when one 
increases the representative algebra dimension, the multiplication result will be closer to the usual Minkowski 
product. We have compared two approaches for interpreting substraction operation, and the canonical 
approach we have proposed, called true arithmetics is more coherent and efficient. Differential calculus 
is possible and very efficient to solve some optimization problems. It is now possible to build inclusion 
functions from the natural ones. This will be studied in a more accurate way in a forthcoming paper. The 
set of intervals is now endowed with an order relation which permits to define inequalities for intervals. One 
of the straightforward application can be non-linear simplex algorithm, the so-called Nelder-Mead simplex 
method or downhill simplex^® ^"^ which is derivative-free method and can be easily implemented. We have 
exhibited some examples of applications : optimization, diagonalization and inversion of matrices which 
clearly state that the arithmetic is stable and that if the initial datas are known with a certain uncertainity 
(belonging to an interval), it is thus possible to estimate with accuracy the point solution of the problem, a 
real number or a smaller interval centered around it. 
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